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Abstract — In previous work, we have developed a dynamic 
learning paradigm for "programming" a general quantum com- 
puter. A learning algorithm is used to find a set of parameters 
for a coupled qubit system such that the system at an initial time 
evolves to a state in which a given measurement results in the 
desired calculation value. This can be thought of as a quantum 
neural network (QNN). Here, we apply our method to a system 
of three qubits, and demonstrate training the quantum computer 
to estimate both pairwise and three-way entanglement. 

Index Terms — quantum algorithm; entanglement; dynamic 
learning 



I. Introduction 



Detection and quantification of entanglement remains an 
important basic challenge in quantum information. Many dif- 
ferent measures have been proposed IH-El. For a system of 
two qubits (quantum bits) there does exist a general measure 
0, but for larger systems [4] there are many problems. Most 
procedures require optimization 0, and/or reconstruction 
of the density matrix (7), both of which can rapidly become 
impractical with increasing size of the system. Also, for 
systems larger than two qubits there exist inequivalent kinds 
of entanglement |8|. The "witness" approach (9) is useful, but 
requires that the state of the system be "close", in some sense, 
to a given, known state. 

Recently Behrman et al. iflOl . ifTTI proposed using an 
approach based on adaptive computing: finding a set of pa- 
rameters for the time evolution of a two-qubit system such that 
the qubit-qubit correlation function at the final time is mapped 
onto an entanglement witness. Remarkably, this witness did 
not require "closeness" to a given initial state. Here, we extend 
our work to a more complex three-qubit system. Our method 
succeeds in finding a single set of parameters such that all three 
qubit-qubit pairwise correlation functions give good witness 
for the pairwise entanglements for large classes of states. 
In addition, the three-point correlation function gives a good 
approximation to the three-way entanglement. No optimization 
procedure or prior state reconstruction is necessary. 
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II. Dynamic Learning for the Coupled Three-qubit 
System: QNN 

A quantum system, whether pure or mixed, evolves in time 
according to the Schrodinger equation: 

d i = -AH,P\ (i) 

at in 

where p is the density matrix and H is the Hamiltonian, whose 
formal solution lfl2l is 

p{t) = exp(iLt)p(0). (2) 

In practice, especially for time-varying Hamiltonians, the 
Schrodinger equation is usually integrated numerically; this 
is the approach we shall use here. We consider a three-qubit 
quantum system whose Hamiltonian is: 

H = K A a xA + K B (t x b + K C (TxC 
+ saOzA + e B o Z B + e c o z c 

+ (aB&zA&zB + (aC&zAVzC + (bC&zB&zC (3) 

where {a} are the Pauli operators corresponding to each of 
the qubits, {K} are the tunneling amplitudes, {e} are the 
biases, and {(}, the qubit-qubit couplings. This represents 
three interacting qubits, labeled A, B, and C, respectively. We 
choose the so-called "charge basis ", which for a system of 
three qubits is the set: 000), |001), |010), |011), |100), |101), 
110), 1 1 1 1 ) . The amplitude for each qubit to tunnel to its 
opposing state, i.e., switch between the and 1 states, is its 
K value; each qubit has an external bias represented by its e 
value; and each qubit is coupled to each of the other qubits, 
with a strength represented by the appropriate £ value. We call 
e a "bias" because, if its value is positive, it will tend to force 
that qubit towards the |0) state (i.e., the energy for being in 
that state is lower), and if its value is negative, towards the |1) 
state. Similarly the sign of ( lowers/raises the energy of the 
corresponding qubits to be aligned/antialigned. Note that, for 
example, a xA acts only on qubit A. We can write that operator 
in the charge basis as an 8x8 matrix: a x ® / <S> I , where a x 

( o 1 \ 

is the familiar 2x2 matrix , / the identity matrix 




, and ® the outer product. Similarly <t z b is the 8x8 



matrix I ® <j z ® I, where a z is the matrix ^ ^ ^ j . Any 

product term, e.g., o z ao z b, can also be written as a single 
matrix: o z ab — <j z ® a z ® I. 

The parameters {K, £,£} control the time evolution of the 
system in the sense that, if one or more of them is changed, the 
way a given state will evolve in time will also change, because 
of Eq. Q] This is the basis for using our quantum system 
as a neural network. There is a mathematical isomorphism 
between Eq. |2] and the equation for information propagation 
in a neural network: 4> output = F W (f> input , where (/) output 
is the output vector of the network, 4>i np ut the input vector, 
and F\y the network operator, which depends on the neuron 
connectivity weight matrix W. Here the role of the "weights" 
of the network is played by the parameters of the Hamiltonian, 
{K, e, (}, all of which can be adjusted experimentally as 
functions of time ff3l . 1141 . By adjusting the parameters 
using a neural network type learning algorithm we can train 
the system to evolve in time to a set of particular final 
states at the final time tf, in response to a corresponding 
set of given inputs. Because the time evolution is quantum 
mechanical (and, we assume, coherent), a quantum mechanical 
function like an entanglement witness can be mapped to an 
observable. Complete details, including a derivation of the 
quantum dynamic learning paradigm using backpropagation 
in time 03], lH6l . are given in ifTTl . 

The time evolution of the quantum system is calculated by 
integrating the Schrodinger equation in MATLAB Simulink 
lTT7l . The ODE4 fixed step size solver was used with a 
integration step size of 0.05 ns. The system was initialized 
(prepared in) each input state, in turn, then allowed to evolve 
for 300 ns. All of the parameters were taken to be functions 
of time; this was done in the simulation by allowing them 
to change to a different constant value every 75 ns (i.e., 
four "time chunks".) Discretization error for the numerical 
integration was checked by redoing the calculations with a 
timestep of a tenth the size; results were not affected. Since the 
error needs to be back propagated through time, the integration 
has to be carried out from tf to 0. To implement this in 
MATLAB Simulink, a change of variable is made by letting 
t' = t f -t. 

III. Pairwise Entanglement Witness 

In earlier work on a two-qubit system, we used a training 
set of only four pure states: the Bell, flat, a correlated product, 
and a partially entangled state. See Table I. Though only 
trained on these four states, the net, when tested on states 
not trained on, successfully reproduced a good approximation 
to the entanglement of formation for large classes of states, 
including mixed states and product states as well as fully and 
partially entangled states. We called this an "experimental" 
entanglement witness, because no prior knowledge of the 
state of the system was required: once the net was trained, 
simple time evolution under the Hamiltonian followed by a 
single measurement was sufficient to detect entanglement for 
any input state. Indeed, the QNN did better than this: under 



most conditions measurement gave an approximate amount of 
entanglement, as well. 

For a three-qubit system, we decided to build on this 
earlier work directly, that is, with an initial training set that 
attempted only to learn the three distinct amounts of pairwise 
entanglement present. Thus, we started with a training set of 
twelve input-output pairs: three copies of the set of four we 
had successfully used earlier for two qubits, one copy for each 
pair of qubits. The input states are shown in Table II, and the 
desired (target) outputs, in Table III. For example, the Bell 
state for qubits B and C (for which we would train the entan- 
glement output for the BC pair to be equal to one) could be any 
(normalized) state of the form: [a |0) + ai|l)] <g) [|00> + |11}]; 
we chose to use the Bell B c state [4||0)]<g)[|00)-|-|ll)]. Table 
IV shows the parameters. The starting values for qubits A and 
B are those found in our earlier work 0T|; these were also 
used as the initial values for qubit C. As before, we took as our 
output the square of the two-qubit correlation function at the 
final time. Thus the target values for the pairwise entanglement 
between qubits B and C were trained to the output function 
(0(t f )) = (a zB (tf)a z c(tf)) 2 = [tr{p{t f )a zB a zC )]\ and 
similarly for the pairs AB and AC. 

Training results are shown in Tables III and IV for the 
outputs and the trained parameters, respectively. RMS error 
per output for the set of 36 outputs (three for each member of 
the training set of 12), after 5000 passes through the training 
set, was 1.215 x 10 -3 . To determine whether the net has 
actually learned, we then tested on a large number of states 
not included in the training set, including fully entangled 
states, partially entangled states, product (unentangled) states, 
and also mixed states. A representative sample is shown in 
Table V. RMS error per output for the set of 39 outputs 
was 3.129 x 10 -2 . While errors are significantly larger, it 
is also clear from the output matrix that there is definite 
separation, in two senses: first, that it is easy to see where the 
pairwise entanglement is (e.g., to distinguish between a state 
with AB entanglement and one with BC entanglement); and, 
second, that it is easy to tell the difference among unentangled, 
partially entangled, and fully entangled states. In other words, 
while this is not, strictly speaking, an entanglement measure, 
it is quite a good witness. 

IV. Three-way entanglement 

Clearly the net has been trained successfully to calcu- 
late an approximation to pairwise entanglement. But can it 
distinguish between pairwise and three-way entanglement? 
To find out, we then tested the pairwise-trained net on the 
GHZ states ^(|000) ± |111>). All the pairwise entangle- 
ments tested as very small: (<r z A{tf)<T Z B(tf)) 2 — 0.0031; 
{ozA{t S )(J z c{tf)) 2 = 0.0032; and (a zB (tf)a zC {tf)) 2 = 
0.0015, while the three-point correlation function at the final 
time was significantly non zero: {<J z A{tf)a Z B(tf)a z c(tf)) 2 = 
0.57. This seemed to indicate that the QNN was able to distin- 
guish between pairwise entanglement, as exhibited in the Bell 
and EPR states, and the inequivalent, three-way entanglement 
of a GHZ state. Perhaps the net could be trained to find yet 



TABLE I 

Training set and data for the two-qubit (AB) system fTTl . The relative amplitudes (for the ket states) are given without 

NORMALIZATION FOR CLARITY. THE BELL STATE IS MAXIMALLY ENTANGLED, AND P IS PARTIALLY ENTANGLED. THE MIDDLE TWO ARE PRODUCT 
STATES (FLAT = (|0) + 1 1)) A (|0) + | 1))b AND CY = |0) A (|0) + 7| l))s ) AND THUS HAVE ZERO ENTANGLEMENT. THE CLASSICAL CORRELATION IS 
COMPUTED AS (a z A (0)<7 zB (0)} . CV IS CLASSICALLY CORRELATED BUT NOT ENTANGLED. THE NEXT THREE COLUMNS SHOW THE ENTANGLEMENT AS 
CALCULATED BY THE METHODS OF BENNETT (JJ AND OF VEDRAL fJ)(USING THE VON NEUMANN METRIC AND THE BURES METRIC.) BOTH DISTANCE 
MEASURES HAVE BEEN NORMALIZED TO UNITY. THE LAST COLUMN SHOWS THE TRAINED VALUES FROM OUR PREVIOUS WORK fTTl . FOR 8 = AND 

7 = 0.5, AFTER 2000 EPOCHS; RMS ERROR = 1.08 X 10~ 5 . 
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0.55 


0.32 


0.46 


0.44317 


0.44317 



TABLE II 

Input matrix for training the three-qubit QNN entanglement witnesses. Each column is an input state, showing amounts of each 

OF THE BASIS STATES. THE FIRST TWELVE STATES ARE THE SET USED FOR THE PAIRWISE WITNESS ; ALL THIRTEEN WERE USED FOR THE SECOND 

TRAINING SET. 
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TABLE III 

QNN ENTANGLEMENT TARGET AND OUTPUT MATRICES FOR TRAINING . EACH COLUMN IS AN INPUT STATE CORRESPONDING TO THE RESPECTIVE 

column in Table II; each row corresponds to a target function for the specified entanglement function: AB, (<t z a(*/)o'zs(*/)) 2 ; 

AC, (f z A(*/)c z c(*/)> 2 ; BC, (crzB(tf)<7 z c(tf)) 2 ; AND ABC, {tTzA{tf)°~zB(tf)o'zc{tf)) 2 - The FIRST THREE ROWS ARE THE TARGETS (DESIRED 
VALUES )FOR THE FIRST TRAINING SET (PAIRWISE ENTANGLEMENT ONLY - THE FIRST TWELVE STATES SHOWN IN TABLE II); THE SECOND THREE, THE 
TRAINED VALUES FOR THAT SET; THE NEXT FOUR, THE TARGETS FOR THE SECOND TRAINING SET (ALL THIRTEEN STATES SHOWN IN TABLE II); THE 
LAST FOUR, THE TRAINED VALUES FOR THAT SET. NONZERO TARGET VALUES ARE BOLDFACED FOR EASY COMPARISON. 
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Trained 


AB 


0.9943 


0.0016 
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0.0003 


0.0018 


0.0001 


0.0007 


0.0001 


0.0002 


0.4399 


0.0008 


0.0003 






AC 


0.0000 
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0.0018 
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0.4385 


0.0016 
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0.4392 
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Trained 


AB 


0.9911 


0.0012 


0.0000 


0.0002 


0.0023 


0.0012 


0.0010 


0.0006 


0.0005 


0.4384 


0.0016 


0.0014 


0.0039 




AC 


0.0002 


0.9964 


0.0000 


0.0015 


0.0001 


0.0018 


0.0000 


0.0017 


0.0001 


0.0049 


0.4386 


0.0021 


0.0001 




BC 


0.0002 


0.0012 


0.9969 


0.0002 


0.0001 


0.0000 


0.0000 


0.0009 


0.0000 


0.0026 


0.0006 


0.4396 


0.0066 




ABC 


0.0012 


0.0000 


0.0028 


0.0001 


0.0014 


0.0002 


0.0000 


0.0001 


0.0003 


0.0006 


0.0024 


0.0007 


0.9883 



another output. Buoyed by hope and encouraged by these 
results, we then proceeded to add one additional output for 
the net to train, the three-point correlation function at the final 
time, ((T Z A(tf)(JzB(tf)<T z c{tf)) 2 , and one additional training 
pair to the input matrix: \GHZ-) = t^(|000) - |111», 
for which we set the target outputs as 0, 0, 0, and 1 (for 
each of the two point functions and the three point function, 
respectively: see Table III.) The net was then retrained, begin- 
ning from the values of the parameters found in the previous 
section, with the new training set of thirteen. Results are 



shown in Table III, and the new values for the parameters 
are shown in Table IV. Testing on product states and on 
pairwise-entangled states was repeated; though errors were 
slightly larger, results were substantially similar. We also 
tested \GHZ + ) = ^ (| 000) + 111)); the outputs for the 
pairwise entanglements and for the three-way entanglement 
were indistinguishable from the outputs for the trained state 
\GHZ-). The net now seems to be able to recognize both 
pairwise and three-way entanglement. 

If we had an experimental realization of the three-qubit 



TABLE IV 

Parameters for entanglement, in MHz. Each column is a different parameter, as labelled. Each parameter is a function of time; 
the first column shows the timeslice number. the first four rows are the initial parameters, taken from our earlier work fill ; 
the next four, the parameters after training (5000 epochs ; rms error = 1.215 x 10 -3 ) to the first training set, of twelve, foriust 

pairwise entanglement (the inputs shown in table ii, using the target values in table iii); the last four, after training to the 
second set, of thirteen (the twelve pairwise states plus one three-way.) training for the second set was done for 5000 additional 

epochs; rms error = 1.801 x 10 -3 .) 





Timechunk 


K A 


K B 


K C 


£A 


£B 


EC 


Cab 


Cac 
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Initial 
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2.3576 


2.3576 


2.3576 


0.10913 


0.10913 


0.10913 


0.04503 


0.04503 


0.04503 




2 


2.3576 


2.3576 


2.3576 


0.03768 


0.06377 


0.06377 


0.10117 


0.10117 


0.10117 
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2.3577 


2.3576 


2.3576 


0.08671 


0.03880 


0.03880 


0.10771 


0.10771 


0.10771 




4 


2.3461 


2.3546 


2.3546 


0.07146 


0.07239 


0.07239 


0.04422 


0.04422 


0.04422 


Set 1 


1 


2.0133 


2.4955 


2.4376 


-0.08676 


-0.44110 


0.32887 


-0.35415 


0.06745 


-0.08465 




2 


2.2238 


2.0328 


2.4116 


-0.11451 


-1.08901 


-0.39170 


0.81660 


0.22532 


-0.27369 




3 


2.4349 


2.2618 


2.4203 


-0.05961 


-0.31644 


-0.04549 


-0.61235 


-0.33585 


0.24552 




4 


2.5836 


2.3423 


2.3853 


0.03040 


0.05885 


-0.02624 


0.28887 


-0.07233 


0.13128 


Set 2 


1 


2.1024 


2.5880 


2.3032 


0.39188 


0.55015 


0.79315 


-0.62404 


-0.26318 


-0.41928 
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2.4013 


2.1874 


2.3051 


-0.04748 


-0.38236 


-0.37284 


0.64560 


0.21584 


-0.09610 
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2.4969 


2.2845 


2.3047 


-0.53340 


-0.48951 


-0.34083 


-0.97289 


-0.92255 


-0.74675 
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2.5908 


2.4842 


2.4072 


-0.31082 


-0.10772 


-0.16904 


0.26100 


-0.17536 


0.01217 



TABLE V 

REPRESENTATIVE QNN PAIRWISE ENTANGLEMENT TESTING : TARGET AND OUTPUT MATRICES. NO ADDITIONAL TRAINING WAS DONE: MERELY, THE 
OUTPUTS WERE EVALUATED USING THE TRAINED PARAMETERS FOR TIME EVOLUTION IN TABLE IV, FOR A NEW SET OF INPUT STATES. EACH COLUMN 
IS AN INPUT STATE. THE FIRST THREE ARE EPR STATES OF THE FORM 4j(|01) ± 1 10)) FOR QUBITS AB, AC, AND BC, TIMES A SUPERPOSITION STATE 

FOR THE THIRD QUBIT: ^=(|0) + |1)); THE NEXT THREE, PRODUCT (UNENTANGLED) STATES (THE FOURTH IS THE STATE 

^7= (| 000) + |001) + |010) + |011) + 1 100) + 1 101) + |110) + I 111)) ; the fifth, the state |000); the sixth, the product state 
0.64(0.8|0) + |1))|1)(|0) + 0.7|1))); the next three, P' states of the form |00) ± 1 1 1) for qubits AB, AC, and BC, times a superposition 

STATE FOR THE THIRD QUBIT: 775 (|0) + |1)); THE NEXT THREE, BELL STATES TIMES A SUPERPOSITION STATE; AND THE LAST, A (COMPLETELY) MIXED 

state: M± = 0.5(|000)(000| ± |111)(111|). Each row corresponds to a value for the specified entanglement function: 

(o" z a(*/)o- 2 b(^/)) 2 , (o" z a(*/)°" z c(*/)) 2 . AND {ozB{tf)ozC(tf)) 2 , RESPECTIVELY. THE FIRST THREE ROWS ARE THE TARGETS (DESIRED VALUES); 
THE LAST THREE, THE CALCULATED OUTPUTS, USING THE PARAMETERS IN TABLE III FOR THE FIRST (PAIRWISE) TRAINING SET. UPPER VS LOWER 
SIGNS MADE ESSENTIALLY NO DIFFERENCE TO THE CALCULATED VALUES. NONZERO TARGET VALUES ARE BOLDFACED FOR EASY COMPARISON. 
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Outputs 


AB 
AC 
BC 


0.9365 

0.0000 
0.0005 


0.0000 
0.7991 

0.0001 


0.0004 
0.0023 
0.9419 


0.0019 
0.0010 
0.0080 


0.0006 
0.0002 
0.0013 


0.0051 
0.0006 
0.0081 


0.4798 

0.0003 
0.0011 


0.0022 
0.4034 

0.0022 


0.0001 
0.0017 
0.3276 


0.9661 

0.0028 
0.0026 


0.0005 
0.9227 

0.0002 


0.0004 
0.0009 
0.7692 


0.0002 
0.0004 
0.0010 



system, we could use these results to find (approximately) a 
pairwise or three-way entanglement for any (initial) state of 
the system, by setting the various parameters to the determined 
values, allowing the system to evolve for the set time, then 
measuring the appropriate correlation function. Of course, any 
real experimental system will have nonzero sources of error, 
like noise and decoherence, which will mean the results of 
the simulation will not be exactly correct. One of the great 
advantages of the neural network approach is its ability to deal 
easily with complications like these: with "online" training, 
the QNN can take these into account in the training phase, 
adjusting the values of the parameters accordingly. 

Within the theoretical simulation, it is a simple matter 
to calculate the entanglement of any state we wish. Some 
interesting results are shown in the two figures. In Figure 1, 
we plot the entanglement, as calculated by the QNN, for the 
state a|000) + /3|001} + |010} + |100) , as a function of both a 
and j3. The graph shows two kinds of pairwise entanglement, 
AB (red) and AC (blue), and three-way entanglement (green). 
The three-way entanglement (green) is zero on the scale of the 
graph for all states plotted; the other pairwise entanglement, 



BC, lies right on top of the AC (blue) entanglement on this 
scale, and so is omitted for clarity. Along the a axis, when 
j3 = 1 (left hand edge), the pairwise entanglements AB (red) 
and AC (blue) are equal (i.e., the surfaces meet); this must 
be true, of course, by symmetry. When a = and (3 = 
(front right corner), this is the EPR state in qubits A and 
B, outerproducted with |0)c; its AB pairwise entanglement 
(red) is maximal while other pairwise entanglements (blue), 
as well as the three-way entanglement (green), are zero. 
When a = and f3 = 1 (front corner), this is the so- 
called "W" state (the three-qubit generalization of the 
EPR state), in which entanglement is shared among all three 
qubits, but in a dissimilar manner to that in the GHZ state 
lfl8ll . The net calculates that the W state has all three kinds 
of pairwise entanglement, in (approximately) equal amounts: 
W zA (t f )a zB (t f )) 2 = 0.3899; (a zA (t f )a zC (tf)) 2 = 0.4129; 
and (<J Z B{tf)<Jzc(tf)) 2 = 0.4042, while the three-way 
(GHZ) entanglement is calculated to be essentially zero: 
(a zA (tf)(T zB (tf)(T z c(tf)) 2 = 1.7628 x 10~ 3 . In the other 
direction, at a = 1 and /3 = (back corner), this is the 
two-qubit partially entangled "P" state for the AB pair, and 



Entanglement of a|000>+ p|001 >+|01 0>+|1 oo> Entanglement Of |000>+ a|1 1 0>+ p|1 1 1 > 




1 010) + 1 100)) , as a function of a and /3, as calculated by the QNN, using |000)), as a function of a and /9, as calculated by the QNN, using the trained 

the trained parameters listed in Table IV. Each color represents a different parameters listed in Table IV. The three-point correlation function, which 

entanglement measure: red, the pairwise entanglement between qubits A and represents three-way entanglement, is in green; pairwise AB entanglement, 

B; blue, that between A and C; and green, three-way entanglement among A, in red. The line along which the AB entanglement is equal to the ABC 

B, and C. entanglement is in black, projected onto the af} plane. 



the net correctly calculates that the AB entanglement (red) is 
about 0.4 while the other pairwise entanglements (blue) are 
approximately zero (see Table III for precise numbers.) 

Figure 2 shows three-way (GHZ) entanglement (green) 
and pairwise AB entanglement (red), for the state a|110) + 
/9 1 111) + 1 000), as a function of both a and (3. For a = and 
(3 = 0, this is the product 000) state, for which pairwise (red) 
and three-way entanglement (green) are both zero; for a = 1 
and f3 = 0, this is the BcIIab state (outer producted with 
|0)c), (|11) + |00))<8> |0), which has maximal AB entanglement 
(red) but zero three-way entanglement (green) ; when a = 
and /8 = 1, this is the GHZ state | 111) + 000), which has zero 
AB entanglement (red) and maximal three-way entanglement 
(green). The two surfaces can be seen to cross each other along 
the black line which is projected onto the a/3 plane; on the 
scale of these calculations it is indistinguishable from a = f3. 

V. Conclusion 

We have shown that a quantum computer of three qubits 
can be "trained" to compute, approximately, its own degree 
of entanglement. Pairwise and three-way entanglement can 
each be determined. No prior state reconstruction or tedious 
optimization procedure is necessary, nor is "closeness" to any 
particular state; rather, any unknown state's entanglement can 
be estimated, whether pure or mixed. 

Further work that extends these results to the four-qubit 
system, with comparisons to the three-tangle [19] and other 
measurements [20|, is in progress. 
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